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Abstract 

Based on the NRQCD factorization formalism, we calculate the inclusive charm production rate 
in T(nS') decay at leading order in the strong coupling constant and the relative velocity v of the 
b quark in the quarkonium rest frame. The branching fractions for T(nS) to charm for n = 1, 2, 
and 3 are all around 7-9%. About 60% of the branching fraction into charm is from annihilation of 
the color-singlet bb pair into 7* cc. Most of the remaining branching fraction is from annihilation 
of the color-singlet bb pair decaying into ccgg. We also compute the momentum distributions of 
the charm quark and charmed hadrons in the decays. The virtual-photon contribution to the 
charm quark momentum distribution is concentrated at the end point while the ccgg contribution 
is broad. The momentum distributions for charmed hadrons are obtained by convolving the charm- 
quark momentum distribution with charm fragmentation functions. This makes the momentum 
distributions for charmed hadrons softer than that for the charm quark. This effect is particularly 
significant in the virtual-photon contribution. 

PACS numbers: 12.38.-t, 12.39. St, 13.20.Gd, 14.40. Gx 
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I. INTRODUCTION 



According to the Nonrelativistic QCD (NRQCD) factorization formalism, an annihilation 
decay rate of a spin-triplet 5'-wave (^<S'i) bottomonium is expressed as an infinite series 
of NRQCD matrix elements with corresponding short-distance coefficients [1] . The NRQCD 
matrix elements, which reflect the long-distance nature of the quarkonium, scale as powers 
of the bottom-quark velocity v in the quarkonium rest frame, which is ~ 0.1. At leading 
order in v, the inclusive decay rate of the T is dominated by the color-singlet spin-triplet 
contribution whose NRQCD matrix element is {Oii^Si))^ — (T|(!?i(^5'i)|T), which is defined 
in Ref . [1] . The subscript 1 on the NRQCD four-quark operator Oi denotes that it is a color- 
singlet operator. Thus, at leading order in the inclusive light-hadronic decay rate of the 
T is expressed in a factorized form: 

r[T ^x\^c, i^iijpi, (1) 

mi 

where X represents all possible light hadronic final states into which T can decay and m;, is 
the bottom-quark mass. The short-distance coefficient Ci, which is insensitive to the long- 
distance nature of the T, can be calculated perturbatively. The dimensions of the matrix 
element (Oi('^iS'i))t is 3 so that Ci is dimensionless. 

At leading order in a^, the dominant color-singlet contribution to Ci comes from 
bbii^Si) — > ggg mode, where the three gluons are attached to the bottom-quark line. Here, 
as is the strong couphng constant and hh-ii^^^^Lj) is the color-singlet hh pair with spin s, 
orbital angular momentum L, and total angular momentum J. The leading contribution 
of order to Ci is known through the orthopositronium decay rate obtained by Caswell, 
Lepage, and Sapirstein [2]. The ordcr-a^ corrections to Ci were calculated by Mackenzie and 
Lepage [3]. This result was confirmed recently by Campbell, Maltoni, and Tramontano [4]. 

In addition to the three-gluon mode, Ci may include the virtual-photon contribution 
from bb-i(^Si) ^ ^ qq-\- X. The decay rate is of order e^e^ci;^, where a is the QED 
coupling constant and is the fractional electric charge of the quark q: = | for an 
up- type quark and — | for a down- type quark. That electromagnetic decay rate may appear 
to be highly suppressed compared to the three-gluon mode of order a^. However, we can 

^ Throughout this paper, we suppress the identifier nS in T(nS'), where n is the radial quantum number, 
unless it is necessary. 
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make a rough estimate of the branching fraction Br[T — > 7* — > qq] by using the measured 
branching fractions for Br[T e+e~] [5]: Br[T ^ 7* ^ qq] NcBt[T e+e~]5^^e^, 
where the sum is over the four flavors of quarks hghter than the bottom and A^c = 3 is the 
number of colors. According to this estimate, Br[T — >• 7* — >• qq] pa 6-8%, which may not be 
negligible. 

At higher orders in v, the NRQCD factorization formula (1) must include additional 
contributions from higher Fock states which involve color-octet pairs bbsC^^^^Lj) as well 
as the color-singlet ones, which are suppressed compared to the leading contribution in 
Eq. (1). The order-v^ and order- relativistic corrections to the color-singlet contributions 
were calculated by Keung and Muzinich [6] and by Bodwin and Petrelli [7], respectively. 
Some of the color-octet contributions were also calculated in Refs. [8-10]. 

Because the T is heavy enough, the decay products may include a pair of charmed 
hadrons. However, unlike the light-hadronic decay mode of the T, there has been little 
previous work on open-charm production in T decay. In 1978, Pritzsch and Streng predicted 
the branching fraction of the decay of T into charm to be a few percents [11], where they 
considered T —>■ ggg* followed by g* —>■ cc. In 1979, Bigi and Nussinov took into account a 
fusion process T — > ccg of order a^, in which a pair of virtual gluons create the cc pair [12]. 
In 1992, ARGUS experiment searched for charm production in direct decays of the T(IS') 
to find only an upper limit of Br'^^''[T(15) ^ L>*(2010)± + X] < 0.019 [13]. 

Recent runs of the CLEO III experiment have produced a large amount of data samples 
at the T(15'), T(25'), and T(35') resonances. The S-factory experiments BABAR and Belle 
have accumulated data for T(25') and T{3S) provided by initial-state radiation. The Belle 
Collaboration has also collected data by running on the T{3S) resonance. With these high- 
luminosity data, one can now indeed study open-charm production in T decay. Very recently, 
some of the authors have calculated the total production rates and momentum distributions 
of the charm quark and charmed hadrons, respectively, in the P-wave bottomonium decays 
XbJ — > c -I- X for J = 0, 1, and 2 by using the NRQCD factorization formalism [14]. 
In order to calculate the momentum distribution of the charmed hadrons they used the 
momentum distribution of charmed hadrons measured by the Belle Collaboration in e'^e~ 
annihilation [15]. 

In this work, as an extension of a previous study [14], we consider inclusive charm pro- 
duction in the spin-triplet S-wsve bottomonium decay. At leading order in v, the dominant 
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mechanism for the decay is a color-singlet channel bbi{^Si) — > ggg* followed by g* — > cc. 
As we have described earlier, the color-singlet mode may have significant virtual-photon 
contribution from bbi(^Si) 7* followed by 7* cc. For inclusive charm production, 
the virtual-photon contribution may have a larger fraction than that in the inclusive light- 
hadronic decay because the rate for T ccgg is suppressed by order ag compared to that 
for T ggg. We consider the virtual-photon contribution as well as the QCD contributions 
from bbi{^Si) ccgg and bbi{^Si) — > ccg^ modes. In the current CLEO III analysis on 
the charmed-hadron (h) momentum distribution in T decay, the virtual-photon contribu- 
tion is subtracted experimentally [16]. We therefore also present the results for the QCD 
contributions after excluding the virtual-photon process. 

This paper is organized as follows. In Sec. II, we present the NRQCD factorization for- 
mulas for the inclusive charm production rate and charm momentum distribution in the 
spin-triplet S-wme bottomonium decay. We also discuss the NRQCD matrix element that 
appear as a long-distance factor in the factorization formula. We calculate the charm-quark 
momentum distribution and the total rate for inclusive charm production in the decay in 
Sec. III. In Sec. IV, in order to provide a theoretical prediction that can be compared with 
CLEO III data, we illustrate the charmed-hadron momentum distributions which are ob- 
tained by convolving the charm-quark momentum distribution with fragmentation functions 
for c — > /i that have been fit to e'^e" annihilation data. Finally, a brief summary of this 
work is given in Sec. V. 

II. CHARM QUARK PRODUCTION IN T DECAY 

In this section, we summarize the NRQCD factorization formula for inclusive charm 
production in T decay. In many aspects there is a large overlap with the formalism for 
calculating the inclusive charm production rate in the P-wave bottomonium decay [14]. In 
this work, we follow the same strategies that were employed in Ref . [14] . For details of the 
formalism, we refer the reader to Refs. [1, 14]. 
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FIG. 1: One of six Feynman diagrams for hb — > ccgg. The other five diagrams are obtained by 
permuting the three gluons attached to the bottom-quark Une. 

A. NRQCD factorization formula 

At leading order in v, the NRQCD factorization formula for the inclusive charm produc- 
tion rate in T decay has an analogous form to that for the light-hadronic decay in Eq. (1): 



where is a dimensionless short-distance coefficient that depends on the mass ratio rric/ rrib 
of the charm quark and the bottom quark and c + X represents all possible states containing 
a charmed hadron, such as , D^, Df, A+, or their excited states. 

At leading order in v and as, the dominant source of C[^^ is the decay of a bbi{^Si) pair 
into ggg*, followed by g* cc. One of six Feynman diagrams of the process is shown in 
Fig. 1 and the remaining five diagrams are obtained by permuting the three gluons attached 
to the bottom-quark line. The decay of a bbi{^Si) pair into 'ygg*, followed by g* cc also 
contributes to C[^\ We call the two processes involving g* — > cc as QCD processes. 

Another similar subprocesses that has the same final state ccgg as the QCD process is 
the decay of bbi{^Si) ggY followed by 7* cc. Because the cc pair created in this 
subprocess is in a color-singlet state, there is no interference between this and the QCD 
process whose cc pair is created in a color-octet state. The contribution of the subprocess 
bbi(^Si) gg'-y* followed by 7* cc to the T decay width is suppressed to the QCD 
contribution by a multiplicative factor 8elela'^N^/a'j./{N^ — 4). For as{'my/2) ^ 0.215 and 
a{mr/2) ^ 1/132, the factor is about 0.09%. Therefore, we neglect this subprocess. 

In addition to the QCD modes, a virtual-photon mode, bbi{^Si) 7*, followed by 
7* — s> cc + X contributes to C['^\ A rough estimate of the branching fraction of the virtual- 
photon channel is about 3%, which is the product of the measured leptonic width r[T — *^ 



c + X] = C- 



1 2 



(2) 
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e+e~], coupling e^, and the color factor Nc- This is comparable to the branching fraction 
Br[T ccgg] predicted by Pritzsch and Streng [11]. In this work, we consider the virtual- 
photon contribution cl^^'^ ^ from bbi{^Si) — 7* ^ cc as well as the QCD contribution 
C[^^^ -* which is composed of bbi{^Si) — >■ ccgg and bbi(^Si) ccg'y. Then the short- distance 
coefficient C['^^ is expressed as 

The two QCD contributions to C['^^^ ^ are essentially the same except for overall factors 
so that 

(ji'^/9*) _ (j{ccgg) _|_ Q{ccgi) _ p q^ccqq) 

■e 

The factor is defined by 



where cf''^ and cf'^^ are contributions of ccgg and ccg'^ channels to c|'^^^ respectively. 



^^-1 + ^^- (^) 

For as{ra-Y/2) a; 0.215 and a{mr/2) ^ 1/132, the numerical value for is about 1.02, 
which indicates that the ccg^i contribution is only about 2% of the QCD contributions. 
Since we are interested in the total decay rate oil ^ c + X and distributions of the charm 
quark with respect to its kinematic variables, we integrate out the phase spaces for the 
gluons or the photon. Therefore, the total and the differential width of the QCD mode are 
the same as those for the ccgg final state up to the multiplicative factor F-y in Eq. (5). 

As in the case of light hadronic decay of the T, inclusive charm production in T decay 
may have contributions from the decay of the color-octet pair fefeg through fofeg ~^ 9* ■, followed 
by g* cc. While the color-octet contribution is suppressed to the color-singlet contribu- 
tion by order v^, the short-distance coefficient of the octet process bbsi^Si) g* ^ cc is 
enhanced by l/a^. Therefore, the color-octet channel may have non- negligible contribu- 
tions, especially for the decay of higher resonances. In this work, we do not consider the 
color-octet contributions. 

In a recent analysis being carried out by the CLEO Collaboration, as a part of experi- 
mental measurement, the virtual-photon contribution is subtracted from the data samples of 
T c+X by scaling the continuum data by an extra factor based on the branching fractions 
Br[T(n5') — > //"^/i"] and i?hadrons but the data include the contribution from bbi{^Si) — > ^ygg* 
followed by g* — > cc [16]. The CLEO III data should, therefore, be directly compared with 
the QCD contributions C[^^^*^ which excludes C{'/^'^ from Cj^^^ in Eq. (3). 
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In the remainder of this paper, we use the same conventions of superscripts to other 
variables as those used in the short-distance coefficients C^^^ \ C^^\ C^^\ and 
For example, T^"'^ = r[T ^ c + X] and T^^^^^) = r[T ^ ccgg]. 



B. Amplitude for hh annihilation into charm 

The short-distance coefficients Cj'\ cJ'/^*\ cf^^\ cf^^\ and Cf^^*^ are calculable 
by using perturbative matching, which involves the computation of the amplitudes for the 
corresponding perturbative short-distance processes such as hhii^Si) — > ccgg. In this section, 
we summarize a way to calculate the annihilation amplitude for the process hhii^Si) — > 
ccgg. Computations of the amplitudes for the remaining processes are analogous to that for 
hhi(^Si) ccgg. 

At leading order in as, the short-distance process for T — c + X is b{p)b{p) 
c{pi)c{p2)g{p3)giP'i) as shown in Fig. 1. The momenta of the b and the b can be expressed 
in terms of the total momentum P and the relative momentum q of the bb pair: 

p ^ iP + q, (6a) 
P ^ kP-q, (6b) 

where the p and the p satisfy the on-shell conditions p"^ — p^ — ml and P ■ q — In the 



rest frame of the bb pair, P — {2Eb, 0) and q — (0, q), where Ei, — \fn^^-^. In general, 
the perturbative amplitude for the 66 annihilation process is expressed as 

v{p)M\i)b\u{p) = ^{M^b\u{p)v{p)), (7) 

where is a matrix that acts on spinors with both Dirac and color indices. The matrix 

Al[66] for the short-distance process 66 — > c{j)\)c(j)-i)gij)-i)g{j)^ is given by 

•^N - 7^^^H^^(Pi)7^"7A^(P2)e?;(P3)er.(P4) 

X b'A(p - P3 - P4)7"A(p - P4)7" ® T'^T^T'^] , (8) 

perm 

where ^pgrm means the summation over the permutations of the three ghions attached to 
the bottom-quark line. is a generator of the fundamental representation for the SU(3) 
color group and a, 6, and c are color indices for the gluons. ei and £2 are polarization vectors 
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for the external gluons with momenta pa and p4, respectively. Note that the expression in 
Eq. (8) is vahd to any order in v. The function A{k) is defined by 

In general, the amphtude (7) contains contributions other than the color-singlet bbi{^Si) 
state, which we want to project out. A convenient way to carry out the projection is to 
replace the spinor product u(p)v(p) in Eq. (7) by the direct product of the color-singlet 
projection operator tti and the spin-triplet projector e^IIg [7, 17, 18], where 

TTi = ^1, (10a) 

^3 = -TT^FT^ -{^ + m,){f+2E,)ri^-m,), (10b) 

where 1 is the SU(3) color unit matrix and e is the polarization four- vector of the bbi{^Si) 
state so that P -e — 0. The projectors (10) are normahzed as Tr[7ri7r|] = 1 and Tr[(e • n3)(e • 
Hs)^] = 4poPo- At leading order in v, the amplitude for the color-singlet spin-triplet 5'-wave 
bb pair can be written as e^Ai[bbi{'^Si)], where 

^^ = Tr(A^[6fe]n^®7ri) . (11) 

q=0 

Because we are working in the leading order in v, we put q — and, therefore, Ef, — nib- 

The amplitude (11) is finite in the soft limits of any external gluons. At higher orders in v, 
infrared divergences arise in this S-wecve amplitude while a P-wave amplitude has an infrared 
divergence at leading order in v. The amplitude (11) is sensitive to the ratio mc/mb. In the 
massless charm-quark limit rric/mb — > 0, the amplitude (11) acquires a coUinear divergence, 
which cancels that arising from the charm-quark loop corrections to the gluon wave function 
for the T — > ggg process [14]. Because the amphtude (11) is free of any infrared and coUinear 
divergences for rric ^ and g = 0, we do not need any regularization scheme and work in 
four space-time dimensions. 

C. Short-distance coefficients 

We proceed to calculate C[^^ in the NRQCD factorization formula (2). The short-distance 
coefficient C[^^ is insensitive to the long-distance nature of the T. Therefore, C[^^ is calculable 
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perturbatively. This can be done by perturbative matching [1, 14]. In order to determine the 
short-distance coefficients we must calculate the annihilation rate for the color-singlet 
spin-triplet hh state by using perturbative QCD. The perturbative analog of the NRQCD 
factorization formula in Eq. (2) for the annihilation rate of the hbii^Si) pair is 

dT[hh{'S,) - c + X] = dC^^ ^^'^'^'\''^^''^\ (12) 

where the factorization formula (12) is written in a differential form. This form is useful for 
our purpose of calculating the momentum distribution of the charm quark. 



1. Calculation of c'f^^ ^ 

The differential annihilation rate of a color-singlet spin-triplet 5"- wave bb state into charm 
through the process hhii^Si) — > ccgg can be expressed as 

dV[hh{'S,) ccgg] = ( i ^ E -^"^H (l^) 

\ ccgg / 

where Ai is the perturbative amplitude (11) for bbi{^Si) — > ccgg, I^'^ is the spin-1 polariza- 
tion tensor for the bbi{^Si), 

pn pu 

I"'' = -r + -p^, (14) 

is the four-body phase space for ccgg, and X^ccgg indicates summation over the spin 
states of ccgg. The factor 1/3 in Eq. (13) comes from averaging over the spin states for the 
U)\i^S\). A factor of 1/2! is multiplied to the four-body phase space because there are two 
identical particles in the ccgg final state, whose phase spaces arc integrated out. 

In order to complete the matching calculation for dC\ , we need to compute the pertur- 
bative NRQCD matrix element: 

(C»i(35i)),5,(350 - 27V,(2E6)2 = + 0{v\ (15) 

Substituting Eqs. (13) and (15) into Eq. (12) and multiplying by in order to take into 
account the ccg^ process as well as ccgg, we find the differential short-distance coefficient 

dC'tl''^: 



'ivE^^-^r)^. (16) 

^ ccgg / 



9 



If we replace in Eq. (16) with unity or — 1, we get the expression for dC^^^^ or 
^(^(ccff7) ^ respectively. 

A parameterization of the four-body phase space (i$4 for the ccgg final state is derived 
in Appendix A: 

2 TT 

where A(a, 6, c) = a^ + ft^ + c^ — 2a6 — 26c — 2ca. In Eq. (17), we have not set v = so that the 
expression can be used for a more general case, xi, ry, and rx in Eq. (17) are dimensionless 
variables defined by 

xi = Ei/Eb, (18a) 
ry = mY/E^, (18b) 
rx = mx/Ef,, (18c) 

where E'l is the energy of the charm quark in the rest frame of the bb pair and the invariant 
masses mx and my are defined by — X'^ — (p2 + Ps + P4)^ and my = = (pg + ^4)^, 
respectively. The ranges of the integration variables are 

Vr'c< xi <1, (19a) 



< ry < V4 - 4:Xi + rc - v^, (19b) 

where 

re = m^/E^ (20) 

(in2 ciOg are solid-angle elements of the charm antiquark with momentum p2 in the 
X-rest frame and the gluon with momentum pa in the y-rest frame, respectively. 

2. Calculation of c[^^"'''^ 

The virtual-photon (7*) contribution C['^^'^ to the short-distance coefficient C['^ in 
Eq. (4) is proportional to that for the leptonic decay of the T{nS): 

dcf^~^ ^ '^ela'5{l - xi)dxi. (21) 
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The short-distance coefficient C^''^ ' is 

dCt'^*^ = elN, (l + ^) X dCt'- 

= |e^e^7Vca'(2 + rc)v^r^5(l-Xi)dxi, (22) 

where the factor ^Jl — Tc is the ratio of the phase space for the cc final state to the massless 
two-body phase space. The QED couphng in Eq. (22) is defined by the running couphng 
constant a{mr) — 1/131, 1/130.9, and 1/130.9 for the radial quantum number n = 1, 2, 
and 3, respectively. 

D. NRQCD matrix elements 

In order to predict the decay widths (2) for the radial quantum numbers n — 1, 2, and 3, 
we must know the numerical value of (Ci)t for each state. In principle, the NRQCD ma- 
trix elements may be calculated in lattice simulations [19, 20]. However, the only available 
NRQCD matrix element for an 5'-wave bottomonium is that for the IS state [20]. Extrapo- 
lating the value in Rcf. [20], which is obtained from an unquenched calculation using two dy- 
namical light quarks, to three light-quark flavors, one obtains (Oi)t(i5) = 3.95 ±0.43 GeV^, 
where we use the same extrapolation method that is given in Ref. [14]. 

The NRQCD matrix element (Oi)t may also be determined by comparing the NRQCD 
factorization formula for r[T — > e^e~] with the empirical values, which are measured with 
uncertainties of order 2%. Recently, a method to resum a class of relativistic corrections to 
5"- wave quarkonium processes has been developed [21, 22], where order-a^ and resummed 
relativistic corrections are included in the NRQCD factorization formula for the leptonic 
decay width. This method has been used to determine the color-singlet NRQCD matrix 
elements for the S-wme charmonium states [22, 23] and to the calculation of exclusive two- 
charmonium production in e'^e~ annihilation [24, 25]. Using the method given in Ref. [22], 
we can determine the NRQCD matrix elements for each nS state: 

(Oi)t(i5) = 3.071°:?^ GeV^, (23a) 
(Oi)t(25) = 1.62ig:iJ GeV^ (23b) 
(0i)t(35) = 1-281H^ GeVl (23c) 

The uncertainties in the matrix elements in Eq. (23) reflect the errors in the lattice string 
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FIG. 2: Distributions of the scaled momentum fraction yi for the charm quark in decays of T(nS) 
forn = 1 (sohd hne), 2 (dashed hne), and 3 (dotted line) by using rrih = 4.6 GeV, as(mx/2) = 0.215, 
0.212, and 0.210, and (C'i)t = 3.07, 1.62, and 1.28 GeV^ for the IS", 25", and 3S states, respectively. 
The distributions in the range yi < 1 are for the QCD contribution dT^'^/^*\ There is a delta 
function at yi = 1 from the virtual-photon contribution dV^'^^'^ ^ in Eq. (22), which is not shown 
in this figure. 

tension, the one-loop pole mass of the bottom quark, the strong coupling constant, and 
the measured values for the leptonic widths as well as corrections of order f ^ that are not 
included in the potential model [22]. The central value of the matrix element (23a) for 
T(IS') is smaller than the lattice estimate by about 22%. In our numerical calculations, we 
use the values for the NRQCD matrix elements given in Eq. (23). 
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III. DECAY RATE 



A. Charm-quark momentum distribution 

In the quarkonium rest frame, the NRQCD factorization formula for the differential 
distribution with respect to the charm-quark energy fraction xi in T decay is easily deduced 
from Eq. (2): 



(24) 



dxi dxi ml ' 

where the differential coefficient is defined by dC[^^ = dC^^^ ^ + dC'f^'^ ^ , which is analogous 
to Eq. (3), and dC'f^^ "* and dC'f^^ ^ are given in Eqs. (16) and (22). Integration over the 
variables dry-, dVL*^^ and dVL\ is implicit in dC'f^^ ^/dxi on the right side of Eq. (24). 

In the T rest frame, the magnitude of the three-momentum of the charm quark is deter- 
mined by the energy fraction xi as Eq. (Allb). Using this relation, we obtain the charm- 
quark momentum distribution from Eq. (24). It is convenient to introduce the scaled mo- 
mentum fraction yi that is defined as the momentum of the charm quark divided by the 
maximum value that is kincmatically allowed in T decay [14]. The physical range of the 
variable is simple: < t/i < 1. By using the relations between xi and y\ 



xi = \/(l - rc)yf + rc, (25a) 



j ,-y>2 fyi 



we obtain the distribution for the scaled momentum fraction yi. 

dr(^) (1 - r^)yi rfr^"^) 



(26) 



dyi ^{1 - r^)yl + dxi ' 

where dT'^'^^dxi is given in Eq. (24). 

In the NRQCD factorization formula (24), the bottom-quark mass appears in common 
for all three S'-wave states. In order to be consistent with a preceding work [14], we use the 
one-loop pole mass mj, = 4.6 GeV for that m(,. The dimensionless short-distance coefficient 
dC[^^/dxi depends on the strong coupling and the ratio Vc- For the strong coupling, we 
use asimr/^) — 0.215, 0.212, and 0.210 for the ratial quantum number 1, 2, and 3 of the T, 
respectively. The ratio Tc depends on md E^. In the nonrelativistic limit ^; — > 0, which we 
are taking in this work, Tc ^ r = [md'^hf' ■ Then the bounds (19) of the variables x\ and 
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TABLE I: Inclusive charm production rate r*^'^) and partial widths T^'^/^*^ and F^'^/'*'*^ in units of 
keV for = 4.6 ± 0.1 GeV, and (Oi)t in Eq- (23). Uncertainties are estimated as stated in 
the text. The partial widths T^^"^^^) and T^'^'^^'^) can be obtained by multiplying T^"^/^*) by factors 
F'^ ^ 0.982 and 1 - F'^ ^ 0.0184, respectively. 



state \ r (keV) 




r(c/7*) 




T(1S) 


1.47 ± 0.36 


2.60 ± 0.65 


4.07 ± 0.75 


T{2S) 


0.83 ± 0.20 


1.38 ± 0.34 


2.21 ± 0.40 


T{SS) 


0.68 ±0.16 


1.09 ± 0.27 


1.77 ±0.32 



ry are determined by r. We take Vc = r = 4m|,/m^, where m/j = 1.87 GeV is the average 
of the masses of the and D~^. For the 5'-wave spin-triplet bottomonium masses, we use 
mr — 9.46 GeV, 10.02 GeV, and 10.36 GeV, for the radial quantum number n — 1, 2, and 
3, respectively [5]. This choice of Tc correctly reflects the physical kinematics [14]. We take 
the numerical values in Eq. (23) for the color-singlet NRQCD matrix elements. 

Our theoretical prediction for the momentum distribution of the charm quark in the 
inclusive T decay is shown in Fig. 2 in terms of dT'^'^^^*^ dyi, where the solid, dashed, and 
dotted lines are the distributions for the 15*, 25*, and 35* states, respectively. Because the 
color-singlet 5'-wave amplitudes for the processes are infrared finite, the distributions in 
Fig. 2 are finite over the whole range of yi. This is different from those for XbJ c + X 
(J = 0, 1, 2), which grow rapidly as yi — > 1 [14]. The yi distributions are broad because of 
the four-body nature of the ccgg (ccgj) final state. The curves have the maximum values 
2.44 keV at yi = 0.51, 1.38 keV at yi = 0.50, and 1.12 keV at yi = 0.50 for T{1S), T{2S), 
and T(35'), respectively. 

As we have stated in Sec. II A, the distributions in Fig. 2 include only bbi{^Si) ccgg 
and bbi{^Si) — > ccg^ contributions to the short-distance coefficient dC[^^ (16). For yi < 1, 
the leading virtual-photon contribution T — > 7* — > cc does not contribute because the 
distribution is proportional to S{l—yi). The sharp peak from the virtual-photon contribution 
dV^^^"'*y dyi at the end point has contributions to the total inclusive charm production rate 
comparable to that of the QCD contributions illustrated in Fig. 2. 
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TABLE II: Branching fractions Br® = r(*)/r[T] for (i) = {c/g*), (c/7*), and (c), where r«'s are 
given in Table I and r[T] is the measured total width [5] of T for radial quantum numbers n =1, 
2, and 3. 



state \ Br (%) 






Br^'^) 


T(15) 


2.72 ±0.67 


4.81 ± 1.21 


7.53 ± 1.39 


T(25) 


2.60 ± 0.67 


4.30 it 1.12 


6.90 ± 1.37 


T(.35) 


3.31 ± 0.87 


5.36 ± 1.41 


8.70 ± 1.74 



B. Total charm production rate 

The inclusive charm production rate in T decay can be calculated by integrating the 
differential rate (24) over xi or the differential rate (26) over yi. In Table I we list the partial 
widths r(^/5*) and F^^/t*) and the total width F^^^. The partial width V^^'a") for the QCD 
process is the sum of F*^'^'^^^^ and F'^^'^'^^'\ whose values can be obtained by multiplying V'^^/^*^ 
by factors fa 0.982 and 1 — 0.0184, respectively. The theoretical uncertainties in 

Table I are estimated as follows. We consider the uncertainties of the pole mass appearing 
in the factorization formula (2) as m;, = 4.6 ±0.1 GeV. As we have discussed earlier, the 
rUc dependence of the dimensionless short-distance coefficients are completely determined 
hy = r = 4:m\y/m\. Wc do not vary r because the measured values for the hadron 
masses do not contribute to errors at the level of accuracy we take into account. We use 
the uncertainties in Eq. (23) for the NRQCD matrix elements. We consider the errors from 
uncalculated order-i>^ and order-cts corrections by multiplying the central values shown in 
Table I by 10% and by as{mr/2) ^ 0.215, respectively. The error bars in the widths 
appearing in Table I are obtained by combining the uncertainties that are listed above in 
quadrature. 

The branching fractions for the decay channels ccgg + ccg'y, 7* cc, and c + X are 
listed in Table II. The numbers are obtained by dividing the numbers in Table I by the 
measured widths r[T(l,5)] = 54.02±1.25 keV, r[T(2,5)] = 31.98±2.63 keV, and r[T(3,5)] = 
20.32 ± 1.85 keV [5]. Prom Tables I and II, we conclude that the virtual-photon contribution 
is actually greater than the QCD contribution for all three S-wave states. 

At the next-to-leading order in a^, the virtual-photon contribution may contribute to the 
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region yi < 1 because of real-gluon emissions. Because the leading-order contribution to 
r{c/7*) is greater than F^"^/^*) by about factors of 1.6-1.8, the order-cts corrections to F^"^/'''*) 
may modify the shape of the yi distributions for yi < 1 by about 30%. 



IV. CHARMED-HADRON MOMENTUM DISTRIBUTION 

In Sec. Ill, we calculated the momentum distribution of the charm quark and the inclusive 
charm production rate in T decay. Since the charm quark hadronizes into one of charmed 
hadrons with a probability of almost 100%, the charm production rate can be interpreted to 
be the sum of the production rates for the charmed hadrons h. The charmed hadrons include 
the D^, , Df, and A+, which are stable under strong and electromagnetic interactions, 
and excited charmed hadrons, whose decay product includes D^, D^, Df, or A+. As 
is discussed in Ref. [14], the momentum distribution of a charmed hadron produced in T 
decay is softer than that of the charm, because of the effect of hadronization. The momentum 
distribution for a charmed hadron h can be obtained by convolving the charm-momentum 
distribution with a fragmentation function for the charm quark to fragment into a h. For 
more details, we refer the reader to Ref. [14] and references therein. 

The fragmentation function Dc^h{z) gives the probability density for a charm quark with 
plus component of light-cone momentum Ei +pi to hadronize into a charmed hadron h with 
light-cone momentum Eh + ph — z{Ei + pi). The fraction z can be expressed in terms 
of scaled light-cone momentum fractions Zi for the charm and Zh for the charmed hadron, 
which are analogous to the scaled momenta yi and yh [14] , where Zi is 

V(l - rc)yf + + yi - rcyi , . 

1 + ^/l-rc 



Then, the fraction z is expressed as 



^ X , (28) 



where the last factor on the right side of Eq. (28) becomes unity if the difference between 
the mass of the charm quark and that of the charmed hadron can be neglected. Within this 
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approximation the momentum distribution of the charmed hadron can be written as [14] 

dV dzh dzi ^ / I ^ dyi 

1~ ^ 1~ Dc-^h{zh zi)— — 

dyh dyh J,^ zi dzi dyi 

r\._.^_ (s/(l-r.)vi + r-^VT^.y.\^ (29) 

^y{l-rc)yl + rc Jy^ \ ^y{l-r,)yj + r, + ^T^i/i J dyi 

where Vc-^hiz) = zDc-^h{z). 

Exphcit parameterizations for Dc-^h{z) for some charmed hadrons h which are currently 
available can be found, for example, in Refs. [15, 26]. Following a previous work on the 
charm production in the P-wave bottomonium decay [14], we refer to the results obtained 
by the Belle Collaboration [15]. They determined the optimal values of the parameters 
for analytic parameterizations of Dc^h{z) for various charmed hadrons by comparing their 
measured momentum distribution in e'^e" annihilation near the center-of- momentum energy 
10.6 GeV with the distributions predicted by Monte Carlo generators and fragmentation 
functions [15]. 

In our numerical analysis, we consider all the charmed hadrons considered in Ref. [15]. 
These cover all the channels being analyzed by the CLEO Collaboration: T(15') h + X, 
where h = D+, D+, D*+, and A+. We use the Kartvehshvili-Likhoded-Petrov (KLP) 
fragmentation function [27] , which was used in the analysis of charmed-hadron momentum 
distribution in Xb decays. The KLP fragmentation function has a simple parameterization 
depending only on the light-cone momentum fraction z: 

D,^h{z)^NhZ^^il-z). (30) 

The optimal values for the parameter etc determined by the Belle Collaboration are etc = 4.6, 
4, 5.6, 5.6, and 3.6 for D+, Df, D*+, and A+, respectively [15]. 

The normalization is determined by the constraint dzDc-^h{z) =Br[c — > h\. Using 
Table X of Ref. [15], one can infer that the inclusive fragmentation probabilities are Br[c 
h] = 0.565, 0.268, 0.092, 0.220, and 0.075 for the D'\ D+, D+, D*+, and A+, respectively As 
a result, we determine the normalization factors for various charmed hadrons as N^jo = 20.9, 
Nd+ ^ 8.04, Njj+ = 4.59, Nd*+ = 11.0, and iV^+ = 1.93. Note that the branching fractions 
for the c — > D° and c include feeddowns from higher resonances D^* and D*^ . 

Substituting Eq. (30) into Eq. (29) and using the parameters hsted above, we evaluate the 
momentum distributions for the charmed hadrons. In Eq. (29), we use Vc — r — i.m\/m\ 
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for Zh where h — Df, and A+ while Tc — r — Am\)/m\ for zi and zd where D — 
and D+. We use = 1.97 GeV, mD*+ = 2.01 GeV, and = 2.29 GeV [5]. 

In Fig. 3, we show the scaled momentum distributions of the charmed hadrons in 
the decays T{nS) — * h + X , where h = (left column) and (right column) for 
radial quantum numbers n = 1, 2, and 3. The distributions for T{nS) ^ h + X, where 
h — Df (left column) and D*~^ (right column) are illustrated in Fig. 4 and those for the 
A+ baryon in Fig. 5. In Figs. 3, 4, and 5, the dotted and dashed curves represent the 
virtual-photon and QCD contributions, respectively. The solid hues are the sums of the 
two contributions. In contrast to the charm-quark momentum distribution, which has a 
virtual-photon contribution only at the end point yi = 1, the virtual-photon contribution 
for a charmed hadron is smeared out to the region y^ < 1. The peaks of the virtual-photon, 
QCD, and total contributions are placed in the ranges 0.78 < yh < 0.84, 0.28 < yh < 0.34, 
and 0.75 < yh < 0.83, respectively. The softening of the QCD contribution can be seen 
by comparing the QCD contribution to the charmed-hadron momentum distributions in 
Figs. 3, 4, and 5 with the charm-quark momentum distributions in Fig. 2. The softening 
of the virtual-photon contribution is even more dramatic, transforming a delta function at 
1/1 = 1 into the virtual-photon contributions in Figs. 3, 4, and 5. Hadronization significantly 
softens the momentum spectrum of charmed hadrons. 

The charmed-hadron production rates in T decay are obtained by integrating over yh in 
Eq. (29) for each hadron. The rates are the same as the area under the sohd lines in Figs. 3, 
4, and 5. The production rates of D+, Df, and A+ in T(1S) decay are 2.25 keV, 
1.06 keV, 0.37 keV, 0.87 keV, and 0.28 keV, respectively, where the contributions of the 
virtual-photon processes amount to about 65-69% of the production rates. In T(2S') decay, 
we get 1.23 keV, 0.58 keV, 0.20 keV, 0.47 keV, and 0.15 keV for each hadron production 
rate, respectively while in T{3S) decay, 0.98 keV, 0.46 keV, 0.16 keV, 0.38 keV, and 0.12 
keV, respectively. For the decays of the 25" state into charmed hadrons the virtual-photon 
contributions add up to about 63-67% while for those of the 35" state to about 63-66%. 
The fractions of the virtual-photon contributions decrease up to about 63% as the radial 
quantum number n of the T{nS) increases. If we compare the momentum distributions of 
charmed hadrons with that of the charm quark, we could observe that hadronization softens 
the distributions. For the QCD contributions, whose charm-quark distribution reaches lower 
end point, the vertical-axis intercepts for the charmed-hadron momentum distributions sig- 



18 



nificantly shift to the positive direction. This results in loss of probabihty. Because the loss 
is minor for the virtual-photon contributions, whose charm-quark distributions are concen- 
trated at = 1, the fractions of the virtual- photon contributions in the charmed-hadron 
production rates are greater than that of the charm quark production rate. As is discussed 
in Ref. [14], this change in normalization is of order r, which is at the level of the error in the 
fragmentation approximation itself, which is derived from QCD by neglecting corrections on 
the order of the square of the quark mass divided by the hard-scattering momentum [28] . 

As we show in Figs. 3, the production rate of D° is more than twice as large as that of 
This reflects the fact that the D° may be produced from the decays of D*^ and D*'^ 
while D'^ has the contribution only from the feeddown from the decay of D*~^ besides direct 
production from the charm quark. As wc show in Fig. 5, the shape of the scaled momentum 
distribution of the A+ is broader than that of the other charmed hadrons. Because A+ 
is the heaviest among the charmed hadrons that we consider, the smearing effects from 
hadronization is most signiflcant in A+ production. 

There are uncertainties in the fragmentation function. Part of the uncertainties can be 
estimated by comparing the results shown above with those obtained by using different pa- 
rameterizations for the fragmentation functions. We also carried out the same calculations by 
using the CoUins-Spiller (CS) fragmentation functions [29] whose parameters are determined 
in Ref. [15] . With the CS fragmentation functions, the heights of the peaks are smaller than 
those from the KLP fragmentation functions by about 8%. This difference is small enough 
to be within our theoretical uncertainties of about 20%. Another important source of the 
uncertainties may come from the Monte Carlo. The estimate of the uncertainties is out of 
the scope of this work. 

V. SUMMARY 

We have studied inclusive charm production in the decay of the spin-triplet bottomonium 
state T{nS) based on the color-singlet mechanism of the NRQCD factorization formalism. 
The branching fractions of the T{nS) into charm are predicted to be (7.5 ± 1.4)%, (6.9 ± 
1.4)%, and (8.7 ± 1.7)% for n = 1, 2, and 3, respectively. The dominant contribution comes 
from the virtual-photon process bbi{^Si) — >• 7* — >• cc, which contributes about 60% of the 
partial width for inclusive charm production in each T{nS) decay. The remaining portion of 
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FIG. 3: Distributions of the scaled momentum y/j in T{nS) h + X for h = (left column) 
and (right column) in units of keV. KLP parameterization is used for the charm fragmentation 
functions [15]. In each figure, solid, dashed, and dotted curves represent the total, QCD (c/g*), 
and virtual-photon (c/7*) contributions, respectively. 
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FIG. 4: Distributions of the scaled momentum in T(nS') h + X for /i = Df (left column) and 
£)*+ (right column) in units of keV. KLP parameterization is used for the charm fragmentation 
functions [15]. In each figure, solid, dashed, and dotted curves represent the total, QCD (c/g*), 
and virtual-photon (c/7*) contributions, respectively. 
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FIG. 5: Distributions of the scaled momentum y/^ in T{nS) + X in units of keV. KLP 

parameterization is used for the charm fragmentation functions [15]. In each figure, sohd, dashed, 
and dotted curves represent the total, QCD (c/g*), and virtual-photon (c/7*) contributions, re- 
spectively. 
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the partial width comes from the QCD contribution, which is dominated by bbi{^Si) — > ggg* 
followed by g* — > cc. The bbi{^Si) — > ccg^ diagrams contribute only about 0.7% of the total 
widths of T{nS) into charm. 

We have also presented the momentum distribution of the charm quark in the T{nS) 
decays. The virtual-photon contribution is localized at the end point {yi = 1). The QCD 
contributions, which have four-body final states, are broad and have peaks near yi — 0.5. By 
convolving the charm momentum distribution with the fragmentation function Dc-^h{z) that 
have been fit to the e^e" annihilation data by the Belle Collaboration, we have computed 
the momentum distribution of charmed hadrons h — D'^, , Df, and produced 
in T(niS') h + X. The resulting momentum distributions of the charmed hadrons are 
significantly softer than that of the charm quark. The large virtual-photon contribution, 
which is localized at the end point {yi — 1) for the charm quark, is smeared out into a 
distribution with a peak near yh ~ 0.8. The order-cts corrections to the virtual-photon 
contributions may enhance the charm-quark momentum distribution in the range yi < 1. 
The cancellation of the infrared divergence at that order may produce a singular distribution 
at the end point, which is similar to that of the P-wave decays. 

In a recent analysis being carried out by the CLEO Collaboration, as a part of measure- 
ment, the virtual- photon contribution is subtracted from the data. In order to facilitate 
comparison with the experimental data, we have provided all the predictions separately 
for the QCD, virtual-photon, and total contributions. Comparison with the experimental 
results will test the leading-order approximations employed in this work. 
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APPENDIX A: THE FOUR-BODY PHASE SPACE 



The four-body phase space o?$4 is defined by 

where Ei and pi are the energy and momentum of the particle i in the final state. For 
T(P) c{pi)c{p2)g{p3)g{p4), pj ^ pI ^ ml, pj ^ pj ^ 0, and = 2Eb. Since we 

are interested in the momentum distribution of the charm quark in the P-rest frame, we 
evaluate d^4 leaving the three-momentum pi unintegrated: 

d^4 - (^^§^^^3(^ ^ P2 +P3+P4), (A2) 

where X — P — pi. The three-body phase space d^s{X p2 +P3+P4) can be expresses as 
a chain of two-body phase spaces: 

d^siX ^P2+P3+P4)^d^2{X ^P2+Y)—^d^2{Y ^P3+P4), (A3) 

where Y — ps + P4 and my is the invariant mass of Y. When a squared amplitude is 
summed over spin states of all the particles in both initial and final states, the squared 
amphtude becomes independent of the solid angle of pi. Integrating over the solid angle of 
Pi, substituting Eq. (A3) into Eq. (A2), and expressing the two-body phase spaces d^2{X — > 
P2 + Y) and d^2{Y p^ + P4) in the X and Y rest frames, respectively, we find that 



= \P\\\P^\^^-i^ dEidmydnidQ* (A4) 
2^^7r'mx 

where mx — {4:E^ — AEbEi+mlY^'^ is the invariant mass of X. The ranges of the integration 
variables Ei and my are given by 

'mc< El < El,, (A5a) 
< my < nix — (A5b) 
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In Eq. (A4), and Ei are the absolute value of the charm-quark momentum and energy 
in the P-rest frame while P2 (pi) and (^^3) ^^e the three-momentum and the sohd- 
angle element of the c (g) in the X (y)-rest frame, respectively. Explicit components of the 
four- vectors pi, and pg are 

pi = (£;i,0,0,|pi|), (A6a) 

P*2 = {E2,\P*2\ sine* COS (f)*,\pl\ sin e*sm(l)l,\p*\ COS 0*), (A6b) 
pI = HpH, \pI\ sine; cos (f)l,\pl\ sine; sin (f)l,\pl\ cose;), (A6c) 

where (6'2, ^2) (^3) ^3) the polar and azimuthal angles of P2 and pi in the X-rest 
and y-rest frame, respectively In Eq. (A6), 

Ipil = {Ef-miy/', (A7b) 
\p*\ = ^\^l\m\,m\,m% (A7c) 

IP3I = ^, (A7d) 

where A(a, 6, c) = -|- 6^ -|- — 2a6 — 26c — 2ca. Note that in Eq. (A6) P2 a^^^d P3 given in 
the X-rest and F-rcst frame while in order to evaluate Eq. (16) it is convenient to express 
them in the P-rest frame. We introduce the boost matrix bJ^i, transforming an arbitrary 
vector k* = (V^, 0) into k = k): 

W = A'^^fc*^ (A8a) 
A% = 1=. (A8b) 



A°i = A'o = (A8c) 



A^ = 5'^ + ^^^, (A8d) 



^ |fc| 

where i, j = 1, 2, 3. The explicit components of p2 ai'c obtained by boosting P2 from the 
X-rest frame to the P-rest frame, where h.^^ is determined by substituting X into k in 
Eq. (A8). To obtain pa, we need two steps. First, we boost P3 from the F-rest frame to 
the X-rest frame where the boost matrix is given by replacing k in Eq. (AS) by Y . The 
components of Y in the X-rest frame are easily obtained by using y = X — p2 valid in any 
frame. Let the obtained four- vector be p^. It is easy to obtain pa by boosting p^ from the 
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X-rest frame to the P-rest frame. Now we can express all the momenta in the P-rest frame 
so that the Lorentz scalars in Eq. (16) can be represented in the P-rest frame. 
It is convenient to introduce dimensionless variables Xi and ry defined by 

xi = EjEb, (A9a) 
ry = my/Eh, (A9b) 

where the ranges of the variables are 

VFc< xi <1, (AlOa) 
< ry < V4 - 4xi + rc - (AlOb) 

rc is the square of the ratio of the charm-quark mass and E^, rc — ml/ E^. The energy and 
momenta £'2, |pi|, IP2I) ^'^'^ \P%\ ^-^e expressed in terms of the variables Xi and ry: 

El = + (Alia) 
2rx 

Ipil = E,{xl-rcY/\ (Allb) 

\Pl\ = ^^HriyY.rc), (Allc) 

\PI\ = Y^6, (Alld) 

where rx — mx/Eb. Substituting Eqs. (A9) and (All) into Eq. (A4), we obtain the four- 
body phase space in terms of Xi and ry: 

= 5/''^'"' ~ ''^^'^'f^'('''^-'y- ''^^ dx,drYdn;dni. (A12) 
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